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Abstract 

The effects of quenched disorder on the overdamped motion of a driven parti- 
cle on a periodic, asymmetric potential is studied. While for the unperturbed 
potential the transport is due to a regular drift, the quenched disorder induces 
a significant additional chaotic "diffusive" motion. The spatio-temporal evo- 
lution of the statistical ensemble is well described by a Gaussian distribution, 

implying a chaotic transport in the presence of quenched disorder. 
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Stochastic models known as thermal ratchets or correlation ratchets |lj , in which a non- 
zero net drift velocity may be obtained from time correlated fluctuations interacting with 
broken symmetry structures ||, have recently received much attention. This interest is due 
to the wide range of possible applications of these models for understanding such systems 
as molecular motors 0, nano-scale friction ||, surface smoothening ||, coupled Josephson 
junctions ||, optical ratchets and directed motion of laser cooled atoms 0, mass separation 
and trapping schemes at the microscale |§. Recently, spatial disorder in thermal ratchets 
has been shown to reduce the characteristic drift speed PJlOl. Little is known, however, 
about the effects of quenched spatial disorder on the regular or diffusive motion in otherwise 
periodic potentials PHT3|- 

Diffusion-like motion is observed in many types of deterministic systems. In particular, 
it has been shown that in deterministic chaotic systems, diffusion can be normal ||14j| , with 



the mean-square displacement (x 2 ) proportional to time t {(x 2 ) ~ t), or it can be anomalous 
I5| , with (x 2 ) ~ t 7 , (enhanced for 7 > 2, dispersive for 1 < 7 < 2), or have a logarithmic 



time dependence (7 = 0) 

In the present work we report on an unusual behavior that occurs in the case of an 
overdamped ratchet subject to an external oscillatory drive: quenched disorder induces 
a normal diffusive transport in addition to the drift due to the external drive. For the 
parameter range considered, this process is observed even for very small perturbations. 
Moreover, this diffusive motion is enhanced by higher values of the quenched disorder. In 
fact, for high enough disorder the diffusive motion is of the same order of magnitude as the 
regular drift. The possibility of having large fluctuations, of the same order of magnitude as 
the average velocity, can be of great importance for a correct interpretation of experimental 
results. This may be of particular importance in studies of friction, in understanding the 
motion of nanoclusters or monolayers sliding on surfaces, as well as for designing particles 
separation techniques. 

We consider the one-dimensional, overdamped motion of a particle (in dimensionless 
units) on a disordered ratchet: 



= cos(x) + fi cos(2a;) + T sin(u?t) + a (1) 

Here, 7 is the damping coefficient, T and u are, respectively, the amplitude and frequency of 
an external oscillatory forcing, and a £ (x) is the force due to the quenched disorder. For the 
present study, £(x) G [—1, 1] are independent, uniformly distributed random variables with 
no spatial correlations, and a > is the amount of quenched disorder. This corresponds to 
a piecewise constant force on the interval [2kir, 2(k + l)vr), k e Z. The unperturbed ratchet 
potential, 

U(x) = — sin(x) — n sin(2x) (2) 
has been the subject of extensive recent studies [|17| -[T9|]. The quenched disorder term [a 7^ 0) 



is expected to give either a more realistic representation of a real substrate or potential 
landscape, or to model fluctuations in DC current amplitude, as for arrays of Josephson 
junctions. 

It is well known jl7|Jl8| that in the absence of quenched disorder (a = 0) there are 
unbounded solutions of Eq. (0), provided that the driving amplitude T is large enough. 
These solutions tend asymptotically to a constant average velocity independent of the initial 



conditions [|TI| . We have identified a set of parameters where, in the absence of disorder, the 
system shows non-zero current (regular transport). Specifically, we have selected 7 = 1.0, 
fi = 0.25, uj = 0.1, and we have studied the behavior for several values of T with T > 1.4 
(see below). 

For a > the periodicity of the unperturbed potential is destroyed as a result of the 
spatial randomness, and solutions of Eq. (fj) begin to show a very complex behavior, in- 
cluding chaotic motion. The chaotic behavior is characterized by the rate of divergence of 
trajectories starting from very close initial conditions, in other words by the leading positive 
Lyapunov exponent. For a given, fixed realization of the quenched disorder, and for several 
values of T > 1.4, we have computed Lyapunov exponents over trajectories starting from 
origin. We have found positive Lyapunov exponents A ranging from 2.55, for T = 1.4, to 3.22 



for T = 1.76, which shows a very strong chaotic behavior. As a consequence of this chaotic 
behavior, the motion of the particle in the perturbed potential should be characterized by 
ensemble averages performed not only over realizations of disorder, but also over the spatial 
distribution of the positions of the particle in a given realization of the quenched disorder. 

Numerical solutions of Eq. ([I]) were obtained using a variable step Runge-Kutta-Fehlberg 
method p0 |. Averages were performed over ensembles of 5000 trajectories starting from 



different initial conditions very close to the origin x = 0. The ensemble described above was 
left to evolve for 1000 external drive periods T, and every 10 periods the positions x(t) were 
stored for further analysis. 

We have first analyzed the motion in a given realization of disorder. In Fig. p] we show 
results for the second moment, C^t) = ((x(t) — (x(t))) 2 ), where (...) means average over 
the ensemble, as a function of the time t, for two different, fixed realizations of quenched 
disorder, and for two values of the disorder parameter a, a = 0.05 (panel (a)), respectively 
a = 0.10 (panel (b)). The most striking feature is the fact that the second moment which 
is zero in the absence of disorder (corresponding to a purely deterministic motion), becomes 
non-zero in the presence of the perturbation. The non-zero second moment confirms the 
chaotic behavior mentioned above, showing a disorder-induced sensitive dependence on the 
initial conditions. It can be seen also that the time-dependence of the second moment is 
very complicated, and it is dependent on both the realization of quenched disorder and the 
amplitude a of disorder. 

In order to perform averages over the realizations of disorder, we have used for each 
trajectory a different random sequence £(x). In this way, the averages over the ensemble 
of trajectories include also averages over realizations of disorder. Figure ^ shows results for 
the first two moments, C\(t) = (x(t)) and C^i) = ({x(t) — (x(t))) 2 ), as a function of the 
time t for T = 1.5 for several values of disorder parameter a. In contrast to averages over 
a given "landscape", in this case both first and second moment show an asymptotic linear 
dependence on time t, C%(t) ~ v(a) t, C 2 (t) ~ D(a) t. We have considered several other 
values 1.40 < T < 1.76, and we have observed the same linear behavior for all a values 



below a threshold value which depends on T. The quenched disorder induces fluctuations in 
the spatial position around the average value in our system, and the dynamics is no longer 
regular, but rather consists in a superposition of regular drift and diffusion-like chaotic 
motion. Moreover, even for reasonably small amounts of disorder, for example a = 0.1, it 
can be seen that these spatial fluctuations are of the same order of magnitude as the first 
moment, so the knowledge of a particular x(t) no longer gives relevant information about the 
position of the center of mass of the distribution, an observation that can be of importance 
in studies of friction, particularly the sliding motion of clusters on surfaces [§J. 

The fluctuations in the position are characterics of chaotic behavior in deterministic 
systems. The description of the initial ensemble is then given by a probability distribution 
function pt{x), whose first two moments are linear in time as we have shown above. We have 
also calculated the higher order cumulants C n (t), for n < 6, and we have found that they 
increase slower than t n l 2 . Therefore, pt(x) is asymptotically a Gaussian, and it is determined 



by the first two moments |18[. In Figure |3] we show the distributions P(z), where z — x— (x), 
for T = 1.50 and two values of disorder parameter a = 0.05 (panel (a)), and a = 0.10 (panel 
(b)), at several times t, and the scaled distributions f(y) = P{z) x yf(t), where y = z/y/i; 
one can see that the distribution is indeed well aproximated by a Gaussian. This asymptotic 
Gaussian behavior also supports the conclusion that the motion is chaotic, as it was shown 
by Jung et al ||18|| . The reason for this chaotic behavior is the existence of discontinuities 



in the velocity at Xk = 2irk, where k is an integer, introduced by the quenched disorder 
perturbation. These random kicks keep into a transitory regime the trajectories that in the 
absence of disorder would have asymptotically converged to the asymptotic constant speed 
state mentioned above. This "mixing" of transitory regimes causes the chaotic behavior, and 
we emphasize again that it is an effect due solely to the perturbation induced by quenched 
spatial disorder. 

For several values of the external drive amplitude T, we have computed from the slopes 
of the first two moments, the drift velocity v(a), and the diffusion coefficient D(a), as 
functions of the amount of disorder a. The results shown in Fig. f| indicate that below a 



(r dependent) threshold value of a the drift is slightly decreasing with increased quenched 
disorder, while the diffusion coefficient is steadily increasing and tends to saturate at high 
amounts of disorder. The fact that disorder has little effect on the drift motion is explained 
by the fact that the drift is a consequence of the initial asymmetry in the potential, and this 
asymmetry is only weakly influenced by small perturbations. We note here that there is no 
decrease in the diffusion coefficient over the range of disorder considered in this study. This 
is in contrast to the decrease of the diffusion coefficient observed in other systems fl9j|l2l. 
The "divergence" of D(a) above the threshold can be understood if we consider the fact that 
v(a) decreases to zero. For large enough a, some of the trajectories in the ensemble become 
bounded, and their contribution to the second moment is proportional to the displacement 
of the center of mass, thus with t 2 . The number of bounded trajectories increases with time, 
as shown by the steady decrease of the drift velocity toward zero. The contribution to the 
second moment (fluctuations) of the t 2 term thus increases in time, and becomes dominant 
at late time, leading to the above mentioned "divergence" of D(a). 

There is a number of experimental situations where small perturbations of a ratchet po- 
tential are relevant, including such systems as surface electromigration ||, dielectrophoretic 
trapping, and particle separation techniques ||. Our preliminary results for the case of 
a non-negligible inertial term in Eq. (fy) are qualitatively similar to the ones for the over- 
damped case, showing disorder induced chaotic diffusion. In this case, however, both the 
"diffusion coefficient" D and the drift velocity v depend on the mass of the particle. Based 
on the similarities mentioned, our results may be relevant for experiments where the mass 
dependence of the drift velocity or diffusion coefficient is essential. The efficiency of a nano- 
scale surface smoothening by an ac field suggested by Derenyi et al |J could be actually 
significantly smaller than theoretically predicted because of the chaotic diffusion, induced by 
the inherent "disorder" of a real surface, superimposed on the net downhill current. On the 
other hand, the rough, imperfect surface of the electrodes in the dielectrophoretic separation 
technique suggested by Gorre-Talini et al || can actually lead to a better efficiency of the 
process by superimposing the chaotic diffusion and drift on top of the thermal, Brownian 
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motion. Moreover, the ac-separation techniques using a two-dimensional sieve discussed by 
Derenyi and Astumian || can be modified in a very natural way to take advantage of the 
inherent imperfection of the two-dimensional structure. This can be done by replacing the 
Brownian diffusion along the drift direction with an additional ac-field along that direction. 
Also, in this way one can fine tune both the drift velocity and the diffusion coefficient along 
the separation direction by a convenient choice of the ac-field parameters. The temperature 
can then be used for an independent tuning of the electrophoretic mobility, thus for the 
transverse displacement. 

In summary we have shown that the addition of small amounts of quenched disorder in 
the equation of motion of a continuous time system induces a strong diffusive motion. In 
addition, we have found that the presence of small amounts of disorder slightly decreases the 
regular current (drift motion), but significantly increases the transport by chaotic diffusion. 
We have shown also that in the presence of disorder the spatial distribution of positions, 
averaged over the realizations of disorder, is described by a time- dependent Gaussian distri- 
bution, which is a signature of chaotic motion. These unexpected results may help in the 
interpretation of experimental results in studies of friction, particularly at the nanoscale, 
as well as in understanding transport processes in molecular motors or designing particle 
separation techniques. 
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FIGURES 

FIG. 1. Second moment (^(i) as a function of time for two different, fixed realizations of 
quenched disorder (solid, respectively dashed lines). The parameters used in Eq. (||) are 7 = 1, 
H = 0.25, T = 1.76, uj = 0.1, and a = 0.05 (panel (a)), respectively a = 0.10 (panel (b)). 

FIG. 2. (a) First moment C\{t) as a function of time for several values of the amount of 
quenched disorder parameter a. From top to bottom, a = 0, 0.05, 0.10, 0.15, 0.20. (b) Second 
moment C2 (t) as a function of time for several values of the amount of quenched disorder parameter 
a. From bottom to top, a = 0, 0.05, 0.10, 0.15, 0.20. The parameters used in Eq. (|l|) are 7 = 1, 
M = 0.25, T = 1.50, and uj = 0.1. 

FIG. 3. Time evolution of the unsealed spatial distributions P(z), where z = x — (x); from top 
to bottom, time t = 20 T, 40 T, . . . , 100 T. The inset shows the scaled distributions P(z) x \ft vs 
zj \ft. The solid line in the insets shows the theoretical, asymptotic Gaussian form. The parameters 
used in Eq. (||) are 7 = 1, /i = 0.25, T = 1.50, u = 0.1, and a = 0.05 (panel (a)), respectively 
a = 0.10 (panel (b)). 

FIG. 4. (a) Diffusion coefficient D{a) as a function of the amount of quenched disorder for 
r = 1.40, 1.50, 1.55, 1.65, 1.76. (b) Drift velocity v(a) as a function of the amount of quenched 
disorder for T = 1.40, 1.50, 1.55, 1.65, 1.76. The parameters are 7 = 1,// = 0.25, u = 0.1. 
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